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D ' Abstract. We study Poisson and operator algebras with the "quasi-linear property" from 

the Heisenberg picture point of view. This means that there exists a set of one-parameter 
' groups yielding an explicit expression of dynamical variables (operators) as functions of 

"time" t. We show that many algebras with nonlinear commutation relations such as the 
Askey- Wilson, g-Dolan-Grady and others satisfy this property. This provides one more 
(explicit Heisenberg evolution) interpretation of the corresponding integrable systems. 
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> ' 1 Classical version 

Assume that there is a classical Poisson manifold with the Poisson brackets (PB) {x,y} defined 
for all dynamical variables x, y belonging to this manifold. Of course, it is assumed that the PB 
. satisfy the standard conditions: 

O " 

OO ! (i) {x, ct\y + a2z} = a\{x, y} + a2{x, z} - linearity (aq, «2 are arbitrary constants); 

(ii) {x,y} = — {y,x} - antisymmetricity; 
(hi) {x, yz} = z{x, y} + y{x, z} - the Leibnitz rule; 
(iv) {{x, y}, z} + {{y, z}, x} + {{z, x}, y} = - the Jacobi identity. 



If one chooses the "Hamiltonian" H (i.e. some dynamical variable belonging to this manifold), 
then we have the standard Hamiltonian dynamics: all variables x(t) become depending on 
an additional time variable t and the equation of motion is defined as 

x(t) = {x(t),H}. (1.1) 

Of course, the Hamiltonian H is independent of t, because H = {H, H} = 0. More generally, 
a dynamical variable Q is called the integral of motion if Q = 0. Clearly, this is equivalent to 
the condition {Q,H} = 0. 

From (II. ip we have x = {{x, H}, H} and more generally 

— = {...{x,H},H},...H} = {x,H^}, 
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where notation {x, H^} means n repeated PB. Using this formula we can write down the formal 
Taylor expansion of an arbitrary dynamical variable in the form 



x(t) =x(0) + t±(0) + ■■■ + * 



n! dt n 

t n 



+ 

t=o 



= x(0)+t{x{0),H} + --- + -{x(0),H^} + --- . (1.2) 

n\ 

Assume now that there are N dynamical variables a;i(0), £2(0), . . . ,xn(0) (not depending on t) 
with the following property: for every k = 1,...,N the PB of the variable x k (0) with the 
Hamiltonian H depends linearly on all Xj(0): 

N 

{x k {0),H} = J2 F UH)x s (0) + $ k (H), k = l,2,...,N, (1.3) 

s=l 

where Fk s (H), <fr k (H) are some functions of the variable H only. Using rules (i)-(iv) for the PB 
we have 

N N 

{{x k (0),H},H} = ^2F ks (H){x s (0),H} = £ F^(H)x s (0) + 

s=l s=l 

where 

N N 

F ( £(H)=Y J F ki {H)F ts {H), ^\H)=J2Fki(H)MH)- 



i=l i=l 



We see that the double PB {{xfc(O), i?}, i?} is again linear with respect to Xj(0). 

In what follows we will use the matrix notation assuming that F(H) is the N x N matrix 
with the entries Fi k {H) and &{H) is the iV-dimensional vector with the components <&i(H). 
Then it is seen that F^ (H) is the square of the matrix F(H) and the vector is obtained 
by the applying of the matrix F(h) to the vector &(H): = F(H)<fr(H) in accordance 

with usual conventions in linear algebra. 

By induction, we obtain the following formula 

N 

{x k (0),H^} = J2F k i : ) (H)x s (0) + <S> k n - 1) (H), (1.4) 

3=1 

where the matrix F^ n \H) means the n-th power of the matrix Fi k and the vector $( n-1 )(ilf) is 
$ (n-l) =F (n-l)$( H y 

Now using (jl.2p we can write down the formula 

t n 



x k (t) = T{x k (0),H^}- 

n=0 

By (|1.4|) we can present the above formula in the form 

N 

x k (t) = ^ F ks (H; t)x s (0) + G k (H; t), (1.5) 
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where the N x N matrix E(H; t) and the iV-dimensional vector G(H, t) are defined as 

E(H;t) = exp(tF(H)), G(H;t) = (^J exp(TF(H))dT^J *(#). (1.6) 
We used the ordinary definition of the exponential function of the matrix: 
exp(tF) = £ ■ 



ni 

n=0 



and 



We see from (jl.5p that every variable x k (t) is a linear function with respect to initial variab- 
les x s (0). 

Proposition 1. Assume that there is an analytical matrix function E(H;t) and an analytical 
vector function G(H;t) such that condition (jl.5p holds for every k = 1,2, . . . , N . Then this 
condition is equivalent to condition (jl.3p . 



Proof of this proposition is elementary. 

So far, we did not concretize the choice of the Hamiltonian H. Now we assume that all 
variables x s (0), s = 1, . . . , N can be chosen as Hamiltonians: if one puts H = x s (0) then we will 
have linearity property (II. 3j) with respect to all other variables Xi, i ^ s. It is naturally to call 
such algebras the quasi-linear Poison algebras. This means that linearity property holds for all 
generators Xi, apart from the Hamiltonian variable x s . 

Quasi-linear algebras possess a remarkable property: if one chooses any basic variable Xj, 
j = 1,2, . . . , N as a Hamiltonian then for time evolution of all other variables Xi, i = 1, . . . , N, 
i ^ j we have linear property 

N 



i(t) = 2_^ £ijs{Xj;t)x s (0) + r]ij(xj;t), i / j, 



s=l 

where the functions £y S (xj;t) and r)ij(xj;t) can be explicitly calculated as was shown above. 

N 

Notation ^ means that we exclude the term with s = j from the sum (note that excluding 

,5 = 1 

value of s coincides with the fixed "number" j of the Hamiltonian). 

In this case we obtain strong restrictions for the expression of the Poisson brackets {xi,Xk}. 
Indeed, from this property we have that for any pair of variables Xi, x^ one has the relation 

N 

{xj,x k } = F iks {x k )x s + $ik(x k ), i,k = 1,2, . . . , N, i^k (1.7) 

8=1 

with some functions Fi ks (xk) and $i k (x k ). These functions cannot be taken arbitrarily. Indeed, 
we can also write down the similar condition but with reversed order of variables Xi, x k : 

N i 

{x k ,Xi} = ^2 F kis (xi)x s + $ki{xi), i, k = 1,2, ...,N, i^k. 
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But {xk,Xi} = — {x{,Xk}, so we have the system of conditions 

N N 



2^ F kis {xi)x s + $ki(xi) + 2^ F iks (x k )x s + <S> ik (x k ) = (li 



which should be valid for any pair x^, x k , i ^ k. 
Moreover, the Jacobi identity 

{{xi,x k },Xj} + {{x k ,Xj},Xi} + {{xj,Xi},x k } = (1.9) 

should be valid for any triple Xi, x k , Xj with distinct values i, j, k. In our case the Jacobi 
identity (jl.9p is reduced to a system of functional-differential equations 

W ijk + W jki + W ktJ = 0, (1.10) 

where 

N ' 

W ijk = {{Xi, X k }, Xj} = & ik (x k ){x k , Xj} + ^ { F iks{x k ){x s , Xj} + F' iks (x k )x s {x k , Xj}) , 

s=l 

where F'{x) means the derivative of the function F{x). 

Conditions (|1.8|) and p.lOp can be considered as a system of nontrivial functional-differential 
equations for unknown functions Fi ks (x k ) and '^(ccfc). 

Of course, there is a trivial solution of these conditions when all functions are constants, i.e. 
F iks {x k ) and &i k (x k ) do not depend on their arguments x k . If, additionally, &i k = then we 
obtain well-known Lie-Poisson algebras with commutation relations 

N 

{X{,X k } — ^ ^ Cj ks X s 
s=l 

with the structure constants Ci ks satisfying standard restrictions following from (jl.8p and (|1.10p 
(see, e.g. [10]). 

Consider the most general finite-dimensional Poisson algebras with nonlinear Poison brackets 

{xi,x k } = h ik (xi, . . . , x N ), (1.11) 

where hi k (x±, . . . ,xn) are smooth functions of N variables x\, . . . , x^. For general theory con- 
nected with these algebras see e.g. [ID]. In order for the variables Xj, x k to form a Poisson 
algebra, the functions hi k should satisfy some strong restrictions. In particular, in [6j necessary 
and sufficient conditions were obtained in case when hi k (x\, . . . , xjv) are quadratic functions. 
We will provide several examples when Poisson algebras of type (jl.lip satisfy conditions (|1.8p 
and (ll.lOp and hence are quasi-linear. General classification of all Poisson algebras with such 
property is an interesting open problem. 

Note that in [10] the so-called semi-linear Poisson algebras (as well as their operator ana- 
logues) were considered. Such algebras are in general nonlinear but they have linearity property 
with respect to some prescribed generators. However the semi-linear algebras introduced in [ID] 
do not possess, in general, the property fll.TH and hence they are not quasi-linear algebras in our 
sense. 
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2 "Quantum" (operator) version 

Now assume that X k , k = 1, 2, . . . are operators which act on some linear space (either finite or 
infinite-dimensional) . 

The time dynamics is defined by the Heisenberg equations: 

B = [H,B], 

where H is an operator called the Hamiltonian and [H, B] = HB — BH is the commutator. 

As in the previous section, assume that there exist N operators X k and a Hamiltonian H 
such that the conditions 

N 

[H, X k ] = Fks(H)X s + <!> k (H), (2.1) 

s=l 

where F ks (H), § k {H) are some functions depending only on the Hamiltonian H. 
Introduce the n-th repeated commutator 

[H, [H,...,[H,B]...] = a d n H B. 

As in the previous section we have 

N 

^ H X k = Y,F£\H)X s + ^- l \H), 

s=l 

where the matrix F^ n \H) means the n-th power of the matrix F^ k and the vector $( n-1 )(.ff) is 
(H) =F( n -V$(H). 

The Heisenberg time evolution is described by the one-parametric group in a standard manner 

t n 

X k {t) = exp(tH)X k exp(-tH) = X k + t[H,X k ] + ••• + - ad^X fc + • • • . (2.2) 

n ! 

Note that in our approach the time variable t can be an arbitrary complex parameter and 
the "Hamiltonian" H need not be a Hermitian operator. In physical applications it is usually 
assumed that H is a Hermitian operator and we need to change t — > it in order to obtain the 
usual Heisenberg picture where the time t is a real parameter. 

Again, as in classical case we obtain that X k {t) is a linear combination of initial operators X s : 

N 

X k (t) = Y J E ks (H;t)X s + G k (H;t), (2.3) 

s=l 

where the N x N operator-valued matrix E(H; t) and the A^-dimensional operator-valued vector 
G(H,t) are defined by the same formulas (jl.6h . 
Note that if the operators X s form a Lie algebra: 

N 

[Xi,x k ] = y^ff| fc Xg, 

8=1 

where g s ik are the structure constants, then any operator Xj = H taken as a Hamiltonian, 
satisfies conditions (|2.ip with F ks the constants (not depending on H) and <fr k = 0, hence we 
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have formula (|2.2p with G k = which is a standard action of the one-parameter Lie group 
corresponding to the generator H = Xy. 

N 

exp(tX i )X fc exp(-tX j ) = Y J E ks {t)X s . 

8=1 

Of course, in the case of the Lie algebra we can construct full Lie group due to linearity property 
(gf k are constants). However, if F ks {H) depend on an operator H then we see that only a set 
of one-parameter groups exists with the linearity property (|2.2|) . Our next problem therefore, 
will be how to construct algebras satisfying the property (|2,ip for different possible choices of 
the operator H. 

The simplest possibility is the same as in the previous section: we demand that any opera- 
tor Xj can be taken as a Hamiltonian H = Xj, j = 1,2, ... ,N. We call such algebras the 
quasi-linear operator algebras. Condition (|2.ip is replaced with a system of conditions 

N 

[X j ,X k } = Y,'Fjks(X j )X s + ^ jk (X J ), k,j = 1,2..., N, k^j. (2.4) 

s=l 

Then for all possible choices of the Hamiltonian H the time evolution (the Heisenberg picture) 
has the linearity property with respect to all operators apart from the Hamiltonian: 

N 

Xj?\t) = exp(tXj)X k expi-tXj) = ' E jks (X j; t)X s + G jk (X r ,t), 

8=1 

where the functions Ej ks (Xj;t) and Gj k (Xj;t) are easily calculated from the functions Fj ks (Xj) 
and & jk (Xj). 

Of course, compatibility analysis of conditions (12, 4h as well as checking of the Jacobi identity 

\[X U X k ],Xj] + [[X^X&Xi] + [[X^X&Xk] =0, i, j, k = 1, 2, . . . , N 

is rather a nontrivial problem even for the case when all operators X k are finite-dimensional. 
Instead, we present concrete examples of the operator algebras possessing quasi-linear property. 
However first we need to generalize our scheme to include the so-called "extension" operators. 

3 Quasi-linear algebras with extension 

Assume that we have a set of iV so-called "basic" operators X±,X2, . ■ ■ X^ and for any opera- 

(i) 

tor Xi, i = 1, 2, . . . , N we have also "a tower" of extension, i.e. a set of Mi > operators 1^ , 
k = 1, 2, ... , Mi, such that we have the commutation relations 

N Mi 

[Xi,X k ] = Y,' F iks {Xi)X s + Y,G iks (Xi)Y^ + $ ik (Xi), 

s=l s=l 

i,k = 1,2, ...,N, i^k (3.1) 

and 

N / Mi 

\ x i, Y k ] \ = E U iks (Xi)X s + ^ Vi ks {Xi)Y® + W ik (Xi), 

8=1 8=1 

i = l,2,...,JV, k = l,2,...,M l (3.2) 
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with some functions Fik s {x), . . . , Wik{x). In what follows we will assume that all these functions 
are polynomials. 

Thus the extension operators Y^f 1 enter the linear combinations in commutators (|3.ip . (|3.2p 

but only operators Xi can be chosen as "Hamiltonians" . In general, commutation relations 

(i) 

between operators Y k are not defined as well as the commutation relations between the Hamil- 
tonian Xi and operators yP' from another "tower" (i.e. when j ^ i). 

From (|3.ip . (|3.2p it follows immediately that repeated commutators with the "Hamilto- 
nians" Xi have the similar structure 

N Mi 

***** = E + E G 2(*)n w + *w w, *, * = i, 2, . . . , jv 

8=1 8=1 

and 

AT Mj 

S=l 8=1 

with polynomials F^(x), . . . , W^ k \x) which can be obtained from Fik s {x), . . . , Wijfc(a;) by ob- 
vious explicit procedures. 

Hence we have an explicit Heisenberg evolution picture: 

N Mi 

e tx >X k e- tx > = ^ 'FiksiXf, t)X s + ]T G tt8 (X<; t)Y« + $ ik (X i; t), i,k = 1,2, ... ,N 

8 = 1 8 = 1 

and 

N f Mi 

e tx iy Wg-tx* = J- c/ ifes (x i; t)X s + ^ F ifcs (X i; t)y W + Wifc^; t), 

8=1 8=1 

where the functions ^ifc s (a;; t), . . . , Wik(x; t) are obtained from the polynomials F^(x),..., 
wjfc\x) by an obvious way, e.g. 

oo n 
n=0 

etc. 

Thus, as in the previous section, we have an explicit linear Heisenberg evolution of the 

(i) 

operators Xi, Y k under the action of N Hamiltonians Xi, i = 1, 2, . . . , N. In contrast to the 

(i) 

previous section, the operators Y k do not in general provide the linear evolution, i.e. they 

cannot be chosen as Hamiltonians. They serve only as an auxiliary tool in our picture. 

(i) 

The total number M of extension operators Y k can be less then Mi + M<i + • • • + Mat, 
because in special cases some of these operators can coincide, as we will see later. 

It will be convenient to call the corresponding algebras the quasi-linear algebras of the type 
(N,M), where ./V is number of the "true Hamiltonians" X\, . . . ,X^ and M is a total num- 

(i) 

ber of extension operators Y k . For quasi-linear algebras without extension we will use the 
symbol (N,0). 

Of course, essentially the same picture is valid in classical case if one replaces the commutators 
[• • • ] with the Poisson brackets {•••}■ 

In the next sections we construct simplest examples of the quasi-linear algebras in both 
classical and quantum (operator) pictures. 
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4 The qr-oscillator algebra 

Consider the simplest case when we have only two operators X, Y (i.e. N = 2) and we demand 
that property (|2.ip will be satisfied when H is chosen as either X or Y, i.e. in this case we have 
the quasi-linear algebra of type (2, 0) without extension. From previous considerations it follows 
that two conditions 

[X,Y]=F 1 (Y)X + $ 1 (Y), [Y,X]=F 2 {X)Y + <1> 2 (X) (4.1) 

should be valid, where Fi(z), &i(z), i = l,2 are some functions. 
In the classical case we will have similar conditions: 

{x,y} =F 1 (y)x + <Z>i(y), {y, x} = F 2 (x)y + $ 3 (x) (4.2) 

with commutators [• • • ] replaced with the Poisson brackets {•■■}• We see that necessary condi- 
tion for compatibility of (|4.ip or (|4.2p is 

F X {Y)X + $i(y) + F 2 (X)Y + $ 2 (X) = 0. (4.3) 

In the classical case (i.e. when X, Y are two independent commuting variables) we have a simple 
functional equation for 4 unknown functions Fi(z), &i(z), i = 1,2. Assume first that the 
function F 2 (x) is not a constant. Then taking two distinct values x = a and x = (5 and 
considering (14. 3D as a system for unknowns F\(y), &i(y) we easily find that these functions 
should be linear in y. Quite similarly, we find that functions F 2 {x), <& 2 (x) should be linear as 
well. After simple calculations we obtain in this case the most general solution 

{x, y} = axy + fax + fay + 7, (4.4) 

where a, fa, fa, 7 are arbitrary constants (a ^ 0). By an appropriate affine transformation of 
the variables x — > £±x + 171, 2/ — *■ €22/ + V2 we can reduce (14. 4j) to a canonical form 

{x,y} = axy - 1 

in the case if the polynomial axy + fax + fay + 7 is irreducible. If this polynomial is reducible 
(i.e. can be presented as a product of two linear polynomials in x, y) then we have the canonical 
form 

{x,y} = axy. 

In quantum case we have correspondingly either algebra of the form 

[X, Y] = aXY + 7, (4.5) 
if 7 7^ 0, or the Weyl operator pair 

XY = qYX. 

Note that relation (|4.5p can be rewritten in other canonical form 

XY - qYX = 1 (4.6) 



with some parameter q ^ 0,1. The pair of operators X, Y satisfying the commutation rela- 
tion (|4.6p is called the ^-oscillator algebra. 
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In this case the Heisenberg evolution becomes very simple. Indeed, consider, e.g. Y as the 
Hamiltonian. We have from (|4.6|) 

[Y, [Y, X}} = u?Y 2 X - uY, 
where uj = 1 — q, and, by induction, 

ad^X = uj n Y n X -uj^Y 11 - 1 , n = l,2,.... 
Hence we have explicitly 

e tY Xe- tY = e" tY X-0(Y;t), 

where 

e uty _ 1 



<j>{y;t) = / e^dr 

Jo uy 

When u — > (i.e. q — ► 1) one obtains 

e tY Xe~ tY =X-t + ujtY{X - t/2) + 0{lu 2 ). 

Here the first term X — t in the right-hand side corresponds to the simple Weyl shift of the 
harmonic oscillator. Indeed, when X, Y satisfy the Heisenberg- Weyl commutation relation 
[X, Y] = 1 then 



jy 



Xe~ yY =X-t 



i.e. the operator X is shifted by a constant under action of the Hamiltonian H = Y . The second 
term cutY(X — t/2) describes a small "(/-deformation" of the Weyl shift. 

5 The Askey— Wilson algebra 

The Askey-Wilson algebra .AW(3) [HJ [28] can be presented in several equivalent forms. One of 
them consists of 3 operators K\, K2, K3 with the commutation relations 

[Kx,K 2 ] = K 3 , (5.1) 
[K 2 ,K 3 ] = 2 P K 2 K l K 2 + ai{tfi, K 2 } + a 2 K 2 + c x K x + dK 2 + 9l , 
[K^Kx] = 2pK 1 K 2 K 1 + a 2 {K u K 2 } + a x K\ + c 2 K 2 + dK x + g 2 

with some constants p, ai, a 2 , c\, c 2 , d, g±, g 2 . The symbol {• • • } stands for the anticommutator: 
{X,Y} = XY + YX. 

The j4VF(3) algebra belongs to quasi-linear type of (2,1). Indeed, we can rewrite it in the 
form 

[K 1 ,K 2 ]=K 3 , (5.2) 
[K 2 , [K 2 ,Ki]] = R 2 (K 2 )K X + Ri(K 2 )K 3 + R (K 2 ), 
[Kx, [K U K 2 }\ = S^Kx)^ + SxiK^Ks + S (Ki), 



where 



R2{x) = —2px 2 — 2a\x — ci, R\{x) = —2px — a±, Rq(x) = —a 2 x 2 — dx — gx, 
S 2 {x) = —2px 2 — 2a 2 x — c 2 , S\(x) = 2px + a 2 , Sq(x) = —axx 2 — dx — g 2 . 
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We see that ^4^(3) algebra in the form (|5.2p has quasi- linear structure with N = 2, M = 1. 
The operators K±, K 2 play the role of the true Hamiltonians whereas the operator K 3 is the 
only extension. 

If p 7^ one can present the same algebra in an equivalent (more symmetric) form in terms 
of operators X, Y, Z [2gl [Ml W\ 

XY-qYX = Z + C 3 , YZ-qZY = X + d, ZX-qXZ = Y + C 2 , (5.3) 

or, in a pure commutator form 

[X,Y] = (q-l)YX + Z + C 3 , [Y,Z] = (q-l)ZY + X + C 1 , 
[Z,X] = (q- 1)XZ + Y + C 2 . 

In this form we have quasi-linear algebra of type (3, 0) without extension, i.e. in this case all 
3 operators X, Y, Z of the algebra play the role of the Hamiltonians. Hence in this case the 
linearity property holds for all operators X, Y, Z. 

It is instructive to derive directly transformation property of, say, operator X under Hamil- 
tonian action of the operator Y. We have 

ad^X = U n (Y)X + V n (Y)Z + W n (Y), (5.4) 

where U n (x), V n (x), W n {x) are polynomials in x. Initial conditions are obvious: 

U (x) = 1, V (x) = W (x) = 0. 

Derive recurrence relations for the polynomials U n (x), V n (x), W n (x). For this goal let as apply 
the operator ady with respect to (|5.4p : 

ad" +1 X = U n (Y)[Y,X] + V n (Y)[Y, Z] = ((1 - q)YU n (Y) + ^KCn) X 
+ ((1 - q~ l )YV n {Y) - U n {Y)) Z - C 3 U n (Y) + q^dV^Y). 

Hence we have 

U n+1 = (1 - q)YU n + q^Vn, V n+1 = -U n + (1 - q~ l )YV n , 

W n+1 = -C 3 U n + C iq ~ l V n . (5.5) 

From recurrence relations (|5.5p it is possible to find polynomials U n , V n , W n explicitly. In 
particular, it is clear that U n (Y) is a polynomial of degree n and V n (Y), W n (Y) are polynomials 
of degree n — 1. 

For the transformed operator X(t) we have 

X(t) = e tY Xe- tY = Ex(t;Y)X + E 2 (t;Y)Z + E (t;Y), 

where 

El(t;50 .f;™ E^-t^ m -±™&l. 

n=0 n=0 n=0 

From (|5.5p we have a system of linear differential equations for unknown functions Ei(t), E 2 {t) 

Ex(t; x) = (l- q)xE 1 + q~ l E 2 , E 2 (t; x) = -E x + (1 - q~ x )xE 2 . 

Initial conditions are E\(0) = 1, E 2 (0) = 0. This system has constant coefficients (not depending 
on t) and hence can be elementary integrated: 

E x {t-x) = an(x)e^* + a 12 (x)e UJ2 ^ t , E 2 (t;x) = a 21 (x)e^ + a 22 (x)^ t , 



Quasi-Linear Algebras and Integrability (the Heisenberg Picture) 



11 



where coi^(x) are roots of the characteristic equation 



J 2 + x(2 - q - q~ 1 )(x - u) + q" 1 = 



and the coefficients a,ik{x) can be found by standard methods. 
The function Eo(t;x) can then be found as 




We see that the Heisenberg evolution is described by elementary functions (linear combinations 
of exponents) in the argument t. But the coefficients in these linear combinations depend on 
operator H. 

Similar expressions are valid for all other possible choices of the Hamiltonian (i.e. H = X or 
H = Z) due to symmetric form of algebra (15. 3ft . 

Thus for p ^ essentially the same ^4VF(3) algebra can be presented in two equivalent forms - 
either as an algebra of type (2, 1) with the only extension or as an algebra of type (3, 0) without 
extension. However if p = then the only type (2, 1) is known (this special case corresponds to 
the so-called quadratic Racah algebra QR(3) (Hj). 

So far, we have no general results in quantum (i.e. operator) case concerning classification 
scheme. In particular, we do not know, whether ^4VF(3) algebra is the only quasi-linear algebra 
of type (2, 1) or (3,0). In the classical case the situation is slightly better: we already showed 
that for the type (2, 0) the only g-oscillator algebra appears. In the next section we consider 
a classification scheme for the type (3,0). 

6 Classical case of type (3,0) 

Consider the classical case of quasi-linear algebras of type (3,0). This means that we have 
3 dynamical variables x, y, z satisfying the Poisson bracket relations 

{y,z} = Fi(x,y,z), {z,x} = F 2 (x,y,z), {x,y} = F 3 (x,y, z), 
where the functions Fi(x,y,z) should be chosen in such way to satisfy the Jacobi identity 

{{x,y},z} + {{y,z},x} + {{z,x},y} = 0. 
It is easily seen that this condition is equivalent to the relation [6] 



where the vector F has Cartesian components (F\, F 2 , F3) and rot is standard differentiation 
curl operator acting on the vector F. 
Relation (16. lj) has an obvious solution 



where V is the gradient operator and Q(x, y, z) is a function of 3 variables. In this case we have 
the defining Poisson relations in the form 



F,rotF) = 0, 



(6.1) 



F = VQ(x,y,z) 



(6.2) 



{y,z} = Q 



{z, x} = Q 



{x,y} = Q 



(6.3) 



where Q x means derivation with respect to x etc. Note that in this case the function Q(x, y, z) is 
the Casimir element of the algebra, i.e. {x, Q} = {y, Q} = {z, Q} = 0. Poisson brackets of such 
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types are sometimes called the Nambu brackets or Mukai-Sklyanin algebras [T7]. However, the 
Nambu bracket (|6.2|) do not exhaust all admissible Poisson algebras with 3 generators. There 
are non-trivial examples corresponding to rot F^O. 

Consider restrictions on the functions Fi(x,y,z) coming from the quasi-linear property. 

On the one side, choosing y to be a Hamiltonian, we have by this property 

{x, y} = F 3 (x, y, z) = $i(y)x + ^(y)z + $ 3 (y)- (6.4) 

On the other side, choosing x to be a Hamiltonian, we have analogously 

{x,y} = ®4(x)y + <S> 5 (x)z + $ 6 (x) (6.5) 

with some functions $i(x), i = 1,2, ... ,6. 

From (|6.4p and (16. 5j) we see that F 3 (x,y,z) should be a polynomial having degree no more 
than 1 with respect to each variable x, y, z and the most general form of this polynomial is 

F 3 {x,y,z) = a 3 xy + (3 31 x + (3 32 y + #33 z + 73. 

Quite analogously, taking functions F±, F2 we obtain 

3 

Fi(x, y, z) = aix k x[ + ^ (3 is x s + 7*, (6.6) 

s=l 

where x\ = x, x 2 = y, x 3 = z and Oj, f3i s , 7, are some constants. As usual in rhs of (16.60 
notation x^x/ means that triple (i, k, I) has no coinciding entries. 

Condition (|6.ip imposes strong restrictions upon the coefficients ai, %■ It is convenient 
to analyze the canonical forms of obtained algebras up to affine transformations x% — > ^x, + r\i 
with some constants £j, 7]% (of course we demand that ^i^2^3 7^ 0). 

(i) If diagonal entries (3u are all nonzero then necessarily ol\ = a 2 = 03 = a. We will 
assume a/0 (otherwise we will obtain the Lie-Poisson algebras). Then the matrix (3 should 
be symmetric: = /3k%- These conditions are also sufficient for validity of the Jacobi identity 
After an appropriate affine transformation we reduce our algebra to the canonical form 

Fi(x, y, z) = ax k xi + x i + j i 

with only 4 independent free parameters: a and 7$, i = 1,2,3. This is exactly the classical 
version of the Askey- Wilson algebra AVK(3) [12j . 

(ii) If one of diagonal entries is zero, say ^33 = 0, then again we have the same necessary and 
sufficient conditions a\ = a 2 = a 3 = a 7^ and (5ik = Phi- The canonical form is 

F± = (XX2X3 + xi + 71, F 2 = ax\x 3 + x 2 + 72, F 3 = ax 2 x\ + 73. 
In both cases (i) and (ii) we have the Nambu-Poisson brackets (|6.3[) with 

Q = axyz + (x 2 + y 2 + z 2 )/2 + 71 x + j 2 y + j 3 z 
for the case (i) and 

Q = axyz + (x 2 + y 2 )/2 + 71X + j 2 y + 73Z 
for the case (ii). 

(hi) Two diagonal entries are zero, say f3 22 = f3 33 = and j5\\ 7^ 0. Then there exist 
2 possibilities: in the first case a\ = a 2 = a 3 = a 7^ and the canonical form of the algebra is 

F± = ax2X 3 + xi + 71, F 2 = axix 3 + 72, F 3 = ax 2 x\ + 73. 
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This is again the algebra of Nambu type with 
Q = axyz + x 2 /2 + j\x + j 2 y + 73^- 

(iv) As in the previous case we have P22 = P33 = and 0u 7^ 0. But now there is a possibility 
of noncoinciding entries Oj. Namely we have 02 = 03 but a\ ^ 02- We will assume that 
«2 a i 7^ 0. Then we have the canonical form 

F\ = 01X2X3 + xi, F 2 = 02X1X3, F 3 = a 2 x 2 x\. 

This case does not belong to the Nambu type. 

(v) There is further degeneration of the previous case. Namely we can allow a\ = 0, a 2 7^ 
or a 2 = 0, a\ ^ 0. In the first case we have the canonical form 

F\ = xi + 71, F 2 = 02x1x3, F 3 = a 2 x 2 xi. 

In the second case the canonical form is 

F 1 = aix 2 x 3 + xi + 71, F 2 = x 3 , F 3 = x 2 . 

(vi) Assume that all diagonal entries are zero (3\\ = (3 22 = (3 33 = 0. Then if a\a 2 a 3 / we 
have the canonical form 

F 1 =aix 2 x 3 , F 2 = a 2 xix 3 , F 3 = a 3 x 2 xi. 

It is interesting to note that obtained algebras of types (i)-(vi) correspond to the list of 
so-called 3-dimensional skew-polynomial algebras introduced by Bell and Smith [22} [TO] . In our 
case we have classical (i.e. Poisson brackets) analogues of the 3-dimensional skew-polynomial 
algebras. 

Concerning other aspects and applications of the classical (Poissonic) version of the AM^(3) 
algebra see e.g. [12]. 

7 Beyond the AW-algebra 

Consider the so-called Dolan-Grady relations [5] for two operators Aq,A\ 



where lo is an arbitrary constant. These relations generate the so-called infinite-dimensional 
Onsager algebra [1] which plays a crucial role in algebraic solution of the Ising model as well as 
of some more general models in statistical physics [18] . 

On the other hand, we can consider the DG-relations from the "quasi-linear" algebras point 
of view. Indeed, introduce the operators 



L4 ,[A),[A),^i]]] 




[^,[^1,^1,^0]]] 




(7.1) 



A 2 = [A ,A!] 



A 3 



Ai 



[A\, [Ai,Aq]]. 



(7.2) 



Then we see that 




2- 



2, 



adj^i = cv 2n A 2 , 
ad^+% = ~co 2n A 2 , 
ad% +1 A 2 = u 2n A 3 , 
ad^ +1 A 2 = -oj 2n A. 



n = 0,l,2,..., 
n = 0,l,2,..., 
n = 0,l,2,..., 
n = 0,1,2,... 
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Prom these formulas we obtain an explicit Heisenberg evolution of the operators Aq, A\, A2 
if the operators Aq or A\ are chosen as Hamiltonians: 



e tA 0Aie -tA =M + sinhM)^ + cothM)-! 



UJ UJ 
3 



e tA 0Me -tA = CQsHu}t)A2 + ^M^ A 

UJ 

e^A e-^ =A - ^M)^ + coth ^ ~ 1 A 4 , 

UJ UJ 1 

e tM Me -t Al = cosHu}t)A2 _ ™h^0^_ 

We see that the operators Aq, A\ can be chosen as Hamiltonians whereas the operators A%, A3, 
A4 are the extension. Hence the DG-relations provide an example of a quasi-linear algebra of 
the type (2, 3). 

Define now the operator 

W = aA + /3Ai + 7A2 

with arbitrary parameters a, (3, 7. It is seen that 

e tAi We -tA! = aAo + pAi + L cosh ( ut ) _ Q SmhM \ M 



cosh(wt) — 1 sinh(u;i) \ 
+ a 7 ^ A 4 . 

UJ Z u J 



Analogously 



ojM We -rA = aAo + 0Ai + /^ 7COsh ( WT ) + ^mh^r) , ^ 



( ^coshfwr) — 1 sinh(o;r)\ , 

+ $ -^r — + 7 — — 4$, 

\ L0 Z UJ J 

where t and r are arbitrary parameters. If one chooses 

q = w7coth(u;i/2), /? = — a;7coth(cjr/2), 
then the terms containing .A3, A4 disappear and we have 

e tA 1We -tA x = aA(j + _ lA2 = e rA QWe -rA ^ 

whence 

TWT~ l = W, 

where 

T = e - TAo e tAl . 

Equivalently, this means that the operator W commutes with the operator T: 
TW = WT. 

As was noted by Davies [1] this commutation relation was crucial in Onsager's solution of the 
Ising model, where the operators T play the role of the transfer matrix. We see that this relation 
follows directly from the quasi-linear property of the DG-relations. 
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There is an obvious generalization of the DG-relations preserving the quasi-linear property. 

Indeed, define again basic operators Aq, A\ and their extensions A 2 , A3, A4 by (|7,2p . But now 
the triple commutators ad^ Ai and ad^ Ao can contain not only operator A2 as in DG-case ()7. 1|) 
but arbitrary linear combinations of the type 

L4 ,[A ,L4 ,Ai]]] = gi(A )A 1 +g 2 (Ao)A 2 +g 3 (Ao)A 3 + go(A ) (7.3) 

and 

[A 1 ,[A 1 ,[A 1 ,A ]]} = f 1 (A 1 )A + f 2 (A 1 )A 2 + f 3 (A 1 )A A + f (A 1 ) (7.4) 

with some polynomials gi{x), fi(x), i = 0, . . . , 3. We will call these relations the generalized DG- 
relations. The ordinary DG-relations corresponds to the choice 50 = fo = 51 = /1 = 53 = /3 = 
and 52 = -f 2 = uj 2 . 

From (|7.3|) and (|7.4|) it is obvious that for any positive integer n the repeated commutator 
ad^ Q Ai has the same structure 

ad^Ai = g^ ) (A )A 1 + g 2 n) (A )A 2 + 9 { r\A )A 3 + (A ) 
and similarly 

a d n Al A = f[ n) (A 1 )A + f ( 2 n \Ar)A 2 + ft ) (A 1 )A 4 + f^(A 1 ) 

with polynomials gf l \x) and (x) can be easily obtained from gi(x) and fi(x). We thus have 
the "Heisenberg solvability" property 

exp(A t)A 1 exp{-A t) = G 1 {A ;t)A 1 + G 2 (A ; t)A 2 + G 3 {A ; t)A 3 + G {A ;t), 

where functions Gi(x; t) can be easily calculated in the same manner as in (|2.3[) . (Similar relation 
holds for the Heisenberg evolution of the operator A±.) 

Consider a special example of these generalized DG-relations. The so-called "tridiagonal 
algebra" proposed by Terwilliger is generated by the two relations [9] 

[A , A 2 Ai + A^l - PAqA^q - j(A A 1 + A X A Q ) - aA x ] = (7.5) 

and 

[A^AlAo + AqA\ - PAxAqAx - 7i(A A 1 + A^Aq) - aiA Q ] = 0, (7.6) 

where (3, 7, 71, a, a\ are some constants. 

The tridiagonal algebra is closely related with the Askey- Wilson algebra AW(3). Indeed, let 
the operators Aq, A\ belong to the AW(3) algebra of the type (|5.ip . Identify Aq = K\, A\ = K 2 , 
A 2 = K 3 and rewrite relations (|5.ip in an equivalent form getting rid of the operator K 3 = A 2 : 



and 



A\A + A A\ + 2(5 - l)AiA Ai + ai(AiAo + A A X ) 

+ a 2 A\ + ciA + dAi + 51 = (7.7) 



AlA Y + AiAq + 2(5 - l^oAiAo + a 2 (AiA + A Q A{) 

+ mAl + c 2 A x + dA + g 2 = 0. (7.. 
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Now applying the operator ad^ to (|7.8p we obtain relation (|7.5p . Analogously, applying the 
operator ad.A 1 to (|7.7|) we obtain relation (|7,6|) . Parameters of the tridiagonal algebra (|7.5|) . (|7.6|) 
are related with parameters of the AVF(3) algebra as follows 

f3 = 2(1 -p), a = -c 2 , a\ = -ex, 7 = -a 2 , 71 = -ax- 

Thus the tridiagonal algebra (|7.5p . (|7.6p follows from the AH^(3) algebra (|5.ip . However the 
reciprocal statement is not valid: the tridiagonal algebra is larger than ^4^(3) [9]. 

It is easily verified that the tridiagonal algebra (|7.5|) . (|7,6|) can be presented in an equivalent 
form: 

[A , [A Q , [A q , Ax]]] = (2 - P){A A 3 - A 2 A 2 ) + 2 7 A A 2 - jA 3 + aA 2 , 

[Ax, [Ax, [Ax,A ]}} = (2 - /9)(Ai^4 + - 2 ll A l A 2 - 7l A 4 - a x A 2 , (7.9) 

where the operators A 2 , A3, A4 are defined as (|7.2p . Comparing (|7.9|) with (|7.3p and (|7.4|) we see 
that the tridiagonal algebra (|7.5p indeed belongs to the quasi-linear type with the coefficients 

9o = 9i = 0, g 2 (x) = ((3 - 2)x 2 + 27X + a, g 3 (x) = (2 - (3)x - 7, 
/o = /i = 0, / 2 (x) = (2-/3)x 2 -2 7l x-ai, / 3 (x) = (2 - 0)x - 71- 

As in the case of the ordinary DG-relations, the tridiagonal algebra is the quasi-linear algebra 
with two basic operators Aq, Ax and 3 extensions A 2 , A 3 , A4, i.e. it has the same type (2,3) as 
the DG-algebra. 

In a special case 7 = 71 = we obtain the so-called q-deformation of the Dolan-Grady 
relations [241125]: 

[A ,[A ,[A ,Ax} q } q -i} =a[A ,Ax], [Ax, [Ax, [Ax, A } q ] q -i] =ax[Ax,A ], (7.10) 

where [X, Y] q = q 1 / 2 XY - q~ x l 2 YX is so-called g-commutator. The parameter q is connected 
with the parameter /3 by the relation /? = q + q~ l - This deformed DG-algebra plays an important 
role in theory of quantum XX Z Heisenberg model, Azbel-Hofstadter model etc [3j [2]. 

It would be interesting to study what is the meaning of the exact solvability of the Heisenberg 
picture for the g-deformed DG-relations (|7.10p in corresponding exactly solvable models. 

8 Quasi-linear algebras and exactly solvable systems 

Recently Odake and Sasaki proposed an interesting approach to exactly solvable classical and 
quantum mechanical models in the Heisenberg picture [151 116] . They noticed that almost all 
exactly solvable one-dimensional quantum models admit exact solution in both Schrodinger and 
Heisenberg pictures. 

Consider, e.g. the one-dimensional model described by the standard one-dimensional Hamil- 
tonian 

H = p 2 /2 + U(x), (8.1) 

where p is the momentum operator and U (x) is a potential. Of course, the standard Heisenberg 
commutation relation between coordinate and momentum operators [x,p] = i is assumed. 

Introduce the operator X = f(x) with some function f(x) and consider the Heisenberg 
evolution of this operator under the Hamiltonian H: 

t 2 (it) 11 
X(t) = exp(iHt)X exp(-iHt) = X + it[H, X] - —[H, [H, X}] + ■ ■ ■ + ad^X + ■■■ . 
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In general for arbitrary Hamiltonian H the structure of the expression ad^X will be too compli- 
cated, and so it is hopeless to find explicit solution in the Heisenberg picture. Assume, however, 
that the Hamiltonian H and the function f{x) satisfy the restriction 

[H, [H, X}} = ad 2 H X = 9l (H)X + g 2 (H)Y + g (H), (8.2) 

with some functions gi(H), i = 0, 1, 2, where we introduce the operator 

Y = [H,X] = -2if'(x)-f"(x). 

Then it is obvious that ansatz (|8.2p (proposed by Odake and Sasaki [13 [16]) leads to an explicit 
solution 

X(t) = t)X + G 2 (H; t)Y + G (H; t), 

where the functions Gi(H;t), i = 0,1,2 have elementary behavior in time t (i.e. they can be 
expressed in terms of trigonometric or hyperbolic functions) . 

Simple considerations (see [15]) lead to the conclusion that g\{H) and go(H) should be linear 
polynomials in H whereas g 2 (H) should be a constant such that 2g 2 = —g[(H) = —ct\ with 
some constant a\. 

It is easy to find (see [H]) that the most general function f(x) satisfying this property is 
a solution of the equation 

f»( x ) = -l(a 1 f(x) + p l ) 

with arbitrary constants a\, (3\. Hence f(x) is an elementary function: it is either a quadratic 
polynomial in x, or a superposition of two exponents (plus a constant). Corresponding po- 
tentials U{x) can also be easily found [15] - they coincide with well-known "exactly-solvable" 
potentials in quantum mechanics: 

(i) singular oscillator: 

U{x) = a\x 2 + a 2 x~ 2 + ao; 

(ii) the Morse potential 

U(x) = aie~ 2x + a 2 e~ x + a ; 

(hi) the Pdschl-Teller potential 

U(x) = a\ sin -2 x + a 2 cos -2 x + ao 

(in the latter case one can replace trigonometric functions with hyperbolic ones that gives hy- 
perbolic analogs of the Poschl-Teller potential) . 

Observe now that two operators H, X together with the third operator Y = [H, X] form 
some quasi-linear algebra. In this case the operator Y plays the role of extension. Indeed, we 
already established the first commutation relation (|8.2p which in our case can be rewritten in 
the form 

[H,Y] = (a 1 H + a )X-^Y + f3 1 H + (3 . (8.3) 
The second commutation relation is verified easily: 



[Y,X] = -f' 2 (x) = ^X 2 + p 1 X + e 



(8.4) 
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with some constant e. Relations (|8,3p and (|8.4p show that we have a quasi-linear algebra of 
type (2, 1) with two Hamiltonians H, X and the only extension Y. This algebra is equivalent 
to so-called quadratic Jacobi algebra Q J(3) which was considered in [8] as a hidden symmetry 
algebra of exactly solvable Schrddinger Hamiltonians (see also |13j). The Jacobi algebra is 
a special case of the algebra j4W(3) [8]. It corresponds to the following choice of the parameters 
p = a.2 = d = in (15. ip . We thus see that the Odake-Sasaki approach [HI [16] for the Schrddinger 
Hamiltonians H (|8.ip leads naturally to the Jacobi algebra QJ(3). 

Consider a difference analogue of the Schrddinger Hamiltonian. We choose 

H = A(s)T + + C{s)T~ + B(s), (8.5) 

where operator H acts on the space of functions f(s) depending on a variable s and operators 
are the standard shift operators 

T ± f(s) = f(s±l). 

The operator H is a second-order difference operator. Hence the eigenvalue equation 
Hf(s) = Xf(s) 

can be considered as a difference analogue of the corresponding one-dimensional Schrddinger 
equation. 

Introduce also the operators X which acts as a multiplication 
Xf(s) = x(s)f(s) 

with some function x(s) which will be called a "grid" and the commutator Y = [H,X]. 
We now would like to find when the relation 

ad 2 H X = [H, Y] = Wi{H)X + W 2 (H)Y + W (H) (8.6) 

holds, where Wi{H), i = 0, 1, 2 are some polynomials in H. If (|8.6p is fulfilled then obviously we 
have an explicit time dynamics of the "grid" operator X under the Heisenberg evolution with 
the Hamiltonian H: 

X(t) = exp(Ht)X exp(-Ht) = Qi(H; t)X + Q 2 (H; t)Y + Q (H; t) 

with functions Qi(H;t) which can be easily explicitly calculated from Wi(H). 

We will assume that A(s)C(s) ^ and the grid x(s) is nondegenerated, i.e. x(s) ^ x(s + 1) 
and x(s) / x(s + 2). Under such conditions it is always possible to assume that C(s) = 1. 
Indeed, using similarity transformation H — ► F(s)HF^ 1 (s) with some functions F(s) one can 
achieve the condition C(s) = 1. Such similarity transformation does not change the operator 
relation (|8.6p . 

Assume that polynomial W%{H) have a degree iV with respect to the variable H. Then 
it is clear from nondegeneracy of x{s) that degrees of polynomials W\{H) and Wq{H) cannot 
exceed N + 1. Assuming that deg(Wi) = deg(Wo) = N + 1 (otherwise we again obtain 
degeneracy of x(s)) we obtain from (|8.6p an operator identity of the form 

N+l 

J2 E k {s)T k = Q 
k=-N-l 

with coefficients Ek(s) which can be explicitly calculated after substitution of (|8.5|) and X = x(s) 
into ()8.6|) . From this identity we should have 



E k (s) = 
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for all k = —N — 1, —N, . . . ,N + 1. Assume that N > 2. Then the highest-order conditions 
-E_tv-i( s ) = Ejv+i(s) = are reduced to very simple first-degree recurrence relations for the 
grid x(s): 

£atx(s + 1) + r] N x(s) + Cat = 

and 

6vx(s - 1) + r/Ar^( s ) + Cn = 0, 

where £at, tjn, are some complex coefficients not depending on s. It is easily verified that 
for all possible choices of the coefficients these conditions are incompatible with nondegeneracy 
conditions x(s) ^ x(s + 1) and x(s — 1) ^ x(s + 1) for the grid x(s). Hence necessarily 
N < 1. Thus W2(H) should be a linear polynomial (or a constant) in H and Wq(H), W\(H) 
should be quadratic (or linear) polynomials in H. Consider again the highest-order conditions 
E 2 (s) = E_ 2 {s) =0: 

£x(s + 1) + r]x(s) + x(s - 1) + C = 

and 

£x(s - 1) + rjx(s) + x(s + 1) + C = 

with some constants £, r/, 

These conditions are compatible with non-degeneracy conditions for the grid x(s) only if 
£ = 1 and then we obtain the linear recurrence relation 

x(s + 1) + x(s - 1) + r)x(s) + C = 0. 

This equation is well known - it describes the so-called ' Askey- Wilson grids" ( AW-grid) [H . 
The most general is the "(/-quadratic" grid: 

x(s) = c\q s + c 2 q~ s + c 

with some constants cq, c±, c 2 and a "basic" parameter q. In the limit q = 1 we obtain the 
quadratic grid: 

x(s) = c 2 s 2 + ClS + Co- 

Further simple analysis of the relation f)8 . 6|) allows to obtain explicit expressions for A(s) 
and B(s). It appears (we omit technical details) that corresponding operator H coincides with 
the difference Askey- Wilson operator on the grid x(s) (for details see [U HU [26]). Thus we 
showed that the only condition (|8.6[) concerning integrability in the Heisenberg picture of the 
second-order difference operator leads uniquely to the Askey-Wilson difference operator and 
corresponding AW-grid. It can be compared with results of [26J, where the similar statement 
was obtained but under different conditions. 

It is then directly verified that the second ("dual") algebraic relation 

[Y,X] = V X (X)H + V 2 (X)Y + V (X) (8.7) 

holds automatically with deg(V 2 (X)) < 1 and deg(V 0) i(X)) < 2. 

Combining (j8.6() and (j8.7[) we immediately conclude that operators H, X together with their 
extension operator Y = [H,X] form the (generic) Askey-Wilson algebra ^4VF(3) (j5.2|) . Thus 
the -AW (3) algebra can be derived uniquely from the "Heisenberg solvability" principle. This 
results seems to be new. Note that in |15| and [16] it was directly verified that the Askey- 
Wilson Hamiltonian (|8.5p (as well as all its special cases corresponding to classical orthogonal 
polynomials) provides an exact solution for the "grid variable" X in the Heisenberg picture. In 
our approach these results follow directly from the quasi- linear property of the AW (3) algebra. 

Concerning classical (i.e. Poisson brackets) analogue of the AW(3) algebra and corresponding 
exactly solvable systems see e.g. [29] . 
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9 Conclusion 

We demonstrated that the "quasi-linear" algebras (with possible extension) can be considered as 
a natural generalization of the Lie algebras: they preserve the linearity property with respect to 
one-parameter subgroups (exponential mapping) constructed from the fixed set of the "Hamil- 
tonians". This provides an exact time evolution (Heisenberg picture) with respect to all these 
Hamiltonians. The situation is almost the same in both classical (with Poisson brackets instead 
of commutators) and quantum picture. Many exactly solvable models in classical and quantum 
mechanics admit an alternative description in terms of a corresponding quasi-linear algebra. 
There are many open questions and problems connected with quasi-linear algebras: 

(i) is it possible to give complete classification of all finite-dimensional quasi-linear algebras 
(with or without extension)? We were able to construct such classification only in the classical 
case (i.e. Poisson algebras) for the simplest cases of types (2, 0) and (3, 0) without extensions. 
We hope that in quantum (operator) case such classification is possible at least in the case when 
all operators are finite-dimensional. 

(ii) in our definition of the quasi-linear algebras, e.g. in condition (12. ip the polynomials F 
(depending on the Hamiltonians H) stand to the left side of the operators X k . Of course it is 
possible to define the "right" quasi-linear algebras with the property 

TV 

[H,X k ] = ^X s F ks (H) + <S> k (H), (9.1) 

s=l 

instead of (12. lj) . It is easily verified that for all considered examples of the quasi- linear algebras 
the "right" version coincides with the "left" version (with possible modification of the structure 
parameters). Is this property valid in general situation? We hope that the answer is positive. 

Note that if the operators H and all Xi are Hermitian and if F ks (x) and Q k {x) are polyno- 
mials in x (as happens e.g. for already considered examples of integrable systems in quantum 
mechanics) then obviously the "left" and the "right" versions are equivalent; indeed, taking 
Hermitian conjugation of (|9.ip we have 

N 

[H,X k ]1 = -[H,X k ] = Y,FL(H)X S + $%(H), 

s=l 

(F* means complex conjugation of corresponding polynomials), i.e. we obtained the "left" version 
from the initial "right" version. However, from general mathematical point of view we cannot 
assume that operators H and Xi are Hermitian, so the problem of equivalence between the "left" 
and the "right" version remains open. 

(hi) possible applications to exactly solvable models. As we know, the Jacobi algebra QJ(3) 
describes exactly solvable one-dimensional quantum Hamiltonians with the singular oscillator, 
Morse and Poschl-Teller potentials. The Askey-Wilson algebra AVF(3) describes all "classical" 
second-order difference equations. By "classical" we mean that these equations have solutions 
in terms of classical orthogonal polynomials on non-uniform grids (the most general are the 
Askey-Wilson difference equations [26]). What about multi-dimensional generalization of these 
results? An interesting approach was proposed in [3], where the q-deformed DG-relations are 
applied to integrable models like XXZ-chwi on the one side and to some generalization of the 
Askey-Wilson polynomials on the other side. 

(iv) the "quasi-linear" algebras do not describe all interesting non-linear algebraic objects in 
theory of integrable systems. Among such algebras are so-called Sklyanin algebras introduced 
in [21] and then generalized in many papers. The Sklyanin algebra consists from 4 generators Si, 
i = 0, . . . , 4 with some special quadratic commutation relations between all possible pairs of the 
generators. These algebras are closely related with elliptic solutions of the Yang-Baxter equation 
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for the 8- vertex model in statistical physics. Recently it was shown that the Sklyanin algebras 
play an important role in theory of biorthogonal rational functions on elliptic grids |2U[ 123] . 
Already on the classical level (Poisson bracket version) the time evolution with respect to any 
"Hamiltonian" Si is described by elliptic functions [llj . Hence in this case time evolution does 
not satisfy linearity property with respect to generators. It would be very interesting to find 
a basic property of these algebras which generalize the quasi-linear property. 
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